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ABSTRACT 

Steady, spherically symmetric, adiabatic accretion and wind flows around 
non-rotating black holes were studied for fully ionized, multi-component fluids, 
which are described by a relativistic equation of state (EoS). We showed that 
the polytropic index depends on the temperature as well as on the composition 
of fluids, so the composition is important to the solutions of the flows. We 
demonstrated that fluids with different composition can produce dramatically 
different solutions, even if they have the same sonic point, or they start with 
the same specific energy or the same temperature. Then, we pointed that the 
Coulomb relaxation times can be longer than the dynamical time in the problem 
considered here, and discussed the implication. 

Subject headings: accretion, accretion disks — black hole physics — hydrody- 
namics — relativity 



1. Introduction 

It is generally inferred from observations tha t the matter fallin g onto black holes is 



of very high temperature, both in microquasars (jCorbel et all 120031 ) as well as in AGNs 



( iRozariska fc Czernyl |2000| ). The electron temperature around 10 9 K and/or the proton 
temperature around 10 12 K or more are accepted as typical values within few tens of the 
Schwarzschild radius, r s , of the central black holes. Moreover, the general theory of relativity 
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demands that the matter crosses the black hole horizon with the speed of light (c) . In other 
words, close to black holes, the matter is relativistic in terms of its bulk speed and/or its 
temperature. On the other hand, at large distances away from black holes, the matter should 
be non-relativistic. 

It is also inferred from o bservations that the astrophysical jets around black hole candi- 
dates have relativistic speeds (IBiretta et a/.ll2003l ). Since the jets originate from the accreting 
matter very close to black holes, their base could be very hot. At a few hundred Schwarzschild 
radii above the disc plane, they can expand to very low temperatures but very high speeds 
(Lorentz factor 7 > a few). And as the fast moving matter of the jets hits the ambient 
medium and drastically slows down to form shocks and hot spots, once again the thermal 
energy increases to relativistic values though the bulk velocity becomes small. Relativis- 
tic flows are inferred fo r gamma-ray bursts (GRBs) too. In the so-called collapsar model 
scenario (jWoosleylll993l ). the collimated bipolar outflows emerge from deep inside collapsars 
and propagate into the interstellar medium, producing GRBs and afterglows. In such model, 
these collimated outflows are supposed to achieve Lorentz factors 7 > 100. 

It is clear in the above examples that as a fluid flows onto a black hole or away from it, 
there are one or more transitions from the non-relativistic regime to the relativistic one or 
vice-versa. It has been shown by quite a few authors that to describe such trans-relativistic 
fluid, the equation of state (EoS) with a fixed adiabatic in dex T {= c p /c v , the ratio of specific 
heats) is inadequate and the relativistically correct EoS rtChandrasekharlll938l ; ISyngdll957l ) 



should be used (e.g., lTaublll948l : Mignone et aZ.ll2005l : iRvu et a/.l 120061 ). 



A fluid is said to be thermally relativistic, if its thermal energy is comparable to or 
greater than its rest mass energy, i.e., if kT > mc 2 . The thermally non-relativistic regime 
is kT <C mc 2 . Here, T is the temperature, k is the Boltzmann constant, and m is the mass of 
the particles that constitute the fluid. So it is not just the temperature that determines a fluid 
to be thermally relativistic, but it is the ratio, T/m, that determines it. Therefore, together 
with the temperature, the composition of the fluid (i.e., either the fluid is composed of 
electron-positron pairs, or electrons and protons, or some other combinations) will determine 
whether the fluid is in the thermally relativistic regime or not. 

The study of relativistic flows around compact objects including black holes was started 
by iMichell (Il972l ). It was basically re casting the t ransonic accretion and wind solutions 
around Newtonian objects obtained by iBondil (jl952j) into the framework of the general the- 
ory of relativity. Since then, a number of authors have addressed t he problem of relativistic 
flows around black holes, each focusing on its various aspects {e.g., 



1976; 


Ferrari 


1985; 


Chakra 


)arti 


1996; 


Das 


2001 


Fukumura & Kazanas 


2007 


Mandal et al. 


2007) 



Blumenthal fc Mathews 



2002; Meliani et al 



2004 



Barai et al. 



20071 ). Barring a few exceptions (e.q.. lBlumenthal &; Mathews 



2006; 



- 3- 



19761 ; iMeliani et all 12004 ) , most of these studies used the EoS with a fixed T, which, as we 
have r ioted, is incapable of describing a fluid from infinity to the horizon. iBlumenthal &: Mathews 
( 119761 ) for the first time calculated the spherical accretion and wind solutions around Schwarzschild 
black holes, wh il e using an ap p roxim ate EoS for the single- com ponent relativistic fluid 
( Mathewd Il97ll ). IMeliani et all (120041 ) modified the EoS used by IBlumenthal fc Mathews 
( 19761 ) to obtain thermally driven spherical winds with relativistic terminal speeds. How- 
ever, there has been no extensive study of the effects of fluid composition on the solutions 
of transonic flows around black holes. We in this paper investigate the effects. 

The paper is organized as follows. In the next section, we present the governing equa- 
tions including the EoS. In section 3, we present the sonic point properties. In section 4, 
we present the accretion and wind solutions. In section 5, we discuss the validity of our 
relativistic EoS. Discussion and concluding remarks are presented in the last section. 



2. Assumptions and Equations 

To ensure that the effects of fluid composition are clearly presented, we keep our model 
of accretion and wind as simple as possible. We consider adiabatic, spherical flows onto 
Schwarzschild black holes. The space time is described by the Schwarzschild metric 

ds* = - ( 1 - cV +(l- 1 dr* + rW + r 2 sin 2 

\ err / \ c A r J 

where r, 6, are the usual spherical coordinates, t is the time, and Mb is the mass of the 
central black hole. Although AGNs and micro-quasars are in general powered by rotating 
flows, studies of spherical flows are not entirely of pedagogic interest. For instance, such 
studies can throw light on the nature of accretions onto isolated black holes in low angular 
momentum and cold clouds. In addition, hot spherical flows may mimic accretions very close 
to black holes, where the accreting matter is expected to be of low angular momentum, hot, 
and with strong advection. Non-conservative processes and magnetic fields are ignored, too. 

The energy-momentum tensor of a relativistic fluid is given by 

= (e + p)v?u v +pg flu , (2) 

where e and p are the energy density and gas pressure respectively, all measured in the local 
frame. The four- velocities are represented by w M . The equations governing fluid dynamics 
are given by 

T£ = and K).„ = 0, (3) 
where n is the particle number density of the fluid measured in the local frame. 
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2.1. EoS for single-component fluids 



Equation (3) is essentially five independent equations, while the number of variables are 
six. This anomaly in fluid dynamics is resolved by a closure relation between e, p and n 
(or the mass density p = nm), and this relation is known as the EoS. The EoS for single- 
component relativistic fluids, which are in thermal equilibrium, has been known for a while, 
and is given by 

e + p K 3 (pc 2 /p) 



pc z 



K 2 (pc 2 /p) 



(4a) 



(jChandrasekharl Il938l ; ISyngd 119571 ) . Here, K 2 and K 3 are the modified Bessel functions of 
the second kind of order two and three, respectively. 

Owing to simplicity, however, the most commonly used EoS has been the one with a 
fixed T, which is written as 

p 



pc 2 + 



r - 1" 



(4b) 



As noted in Introduction, this EoS, whi ch admits the super l uminal sound sp eed, is not 
applicable to all ranges of temperature (IMignone et all 120051 ; iRyu et al\ 120061 ). Here, we 
adopt an approximate EoS 

'9p + 3pc 2 ' 



pc +p 



(4c) 



3p + 2pc 2 / 

which reproduces ve ry closely the re lativistically correct EoS in equation (4a), better than 
the one proposed by iMathewd (119711 ) 



P 



pc 
~3~ 



pc A 



pc_ 

e 



(4d) 



A comparative stu dy of various EoS's for single- component relativistic fluids was presented 



in 



Rvu et all fl2006h . 



2.2. EoS for multi-component fluids 

We consider fluids which are composed of electrons, positrons, and protons. Then the 
number density is given by 

n = Snj = n e - + n e + + n p +, (5a) 

where n e - , n e + , and n p + are the electron, positron, and proton number densities, respectively. 
Charge neutrality demands that 

n e - = n e + + n p + ^> n — 2n e - and n e + = n e -(l — ^), (5b) 
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where £ = n p +/n e - is the relative proportion of protons. The mass density is given by 

1 

7/ 



p = Uriirrii = n e -m e < 2 - M 1 



(5c) 



where 77 = m e /m p , and m e and m p are the electron and proton masses, respectively. For 
single-temperature fluids, the isotropic pressure is given by 



p = Y^pi = 2n e - kT. 



(5d) 



As our EoS for multi-component fluids, we adopt 

'9pi + SriimiC 2 ^ 



e = Se } - = £ 



11^0 +pi 



3pi + 2nj7«jC s 



(5e) 



The non-dimensional temperature is defined with respect to the electron rest mass energy, 
= kT / (m e c 2 ). With equations (5a) - (5d), the expression of the energy density in equation 
(5e) simplifies to 

e = n e -m e c 2 f, (5f) 

where 

/90 + 3/V 



/ = (2-0 



1 + 



90 + 3 
30 + 2 



77 +9 V30 + 2/r7 



The expression of the polytropic index for single-temperature fluids is given by 



where 



p dT 2d& 

e t 2 urr 9kT + 3m lC 2 

®i = — = + kT—— 

ru 3kT + 2m iC 2 



(5g) 



(5h) 



(5i) 



is the energy density per particle of each component. The effective adiabatic index is calcu- 
lated by 



r = 1 + 



N' 



The definition of the sound speed, a, is 



Tp 



2TQ 



c 2 e + p f + 2Q' 



(5j) 



(5k) 



The polytropic index iV (and also the adiabatic index T) is an indicator of the thermal 
state of a fluid. If N — > 3/2 (or T — > 5/3), the fluid is called thermally non-relativistic. 
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On the other hand, if N —>■ 3 (or T — > 4/3), it is called thermally relati vistic. For single - 



component fluids, N and T are given as a function of the temperature alone (IRyu et a/.ll2006l ). 
For mult i- component fluids, however, not just the temperature, the mass of the constituent 
particles also determines the thermal state. Hence, the proton proportion, £, enters as a 
parameter too. 

In Figure 1, we show various thermodynamic quantities and their inter-relations for 
fluids with different £. In Figure la which plots iV as a function of T, the left most (solid) 
curve represents the electron-positron pair fluid (£ = 0) (hereafter, the e~ — e + fluid) and the 
right most (dotted) curve represents the electron-proton fluid (£ = 1) (hereafter, the e~ — p + 
fluid). In the e~ — e + fluid, N — > 3 for kT > m e c 2 , while in the e~ — p + fluid, N — > 3 for 
kT > m p c 2 . In the intermediate temperature range, m e c 2 < kT < m p c 2 , N decreases (i.e., 
the fluid becomes less relativistic) with the increase of £. It is because if £ increases (i.e., 
the proton proportion increases), the thermal energy required to be in the relativistic regime 
also increases. By the same reason, at the same T, the local sound speed, a, decreases as £ 
increases, as shown in Figure lb. However, in Figure lc, it is shown that the relation between 
iV and a is not as simple as the relation between N and T. At the same a, N is smallest for 
the e~ — e + fluid, and it increases and then decreases as £ increases. The behavior can be 
understood as follows. At the same a, as £ increases, the thermal energy increases, but at 
the same time, the rest mass energy increases as well. As noted in Introduction, it is not the 
thermal energy, but the competition between the thermal energy and the rest mass energy 
that makes a fluid relativistic. Consequently, for most values of a, iV increases for £ < 0.2 
and then decreases for £ > 0.2. For very low a, N increases up to £ ~ 0.5, and for very 
high a, N increases up to £ < 0.1. In summary, at a given temperature, the e~ — e + fluid is 
most relativistic, but at a given sound speed, the e~ — e + fluid is least relativistic and fluids 
with finite proton proportions are more relativistic. 



2.3. Equations of motion 

The energy-momentum conservation equation [the first of equation (3)] can be reduced 
to the relativistic Euler equation and the entropy equation. Under the steady state and 
radial flow assumptions, the equations of motion are given by 

r du r 1 / 2 _ _\ 1 dp ,„ . 

U ~JT + 3 = - 1 - Z + u u ( 6a ) 



dr r 2 \ r J e + p dr 

and 



de e + pdn 
dr n dr 
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along with the continuity equation [the second of equation (3)] 



1 dn 
n dr 



2 

r 



1 du r 
u r dr 



(6c) 



Here, we use the system of units where G = Mb = c = 1, so that the units of length and 
time are r g = GMb/c 2 and t g = GMb/c 3 . It is to be noted that in this system of units, 
the Schwarzschild radius or the radius of the event horizon is r s 
calculations, equations (6a) - (6c) are then simplified to 

dv _ (1 -v 2 )[a 2 (2r-3) - 1] 
dr 



r(r — 2)(v — a 2 /v) 



2. After some lengthy 



(7a) 



and 



dQ 
dr 





N 



2r - 3 



+ 



dv 



r(r — 2) v (1 — v 2 ) dr 
where the radial three- velocity is defined as v 2 = —u r u r /(u t u t ). 



(7b) 



For flows continuous along streamlines, equations (7a) - (7b) admit the so-called reg- 
ularit y condition, or the critical point condition, or the sonic point condition (jChakrabarti 



19901 ) that is given by 



and 



1 



la) 



(8b) 



c 2r c - 3 

Here, r c is the sonic point location. Hereafter, the quantities with subscript c denote those 
at r c . From equation (5k), w e know a maiX = l/\/3 (also se e Figure lb). Therefore, from 
equation (8b), we have r c > 3 ( IBlumenthal fc Mathewslll976l ). Since dv/dr = M /V — > 0/0 
at r c , (dv/dr) rc is obtained by the l'Hospital rule 



dv 
dr 



(dU/dr) r 
(dV/dr) r , 



(8c) 



where M and T> are the numerator and denominator of equation (7a). The above equation 
simplifies to 

dv" 



A 



where 



A 



dr 



N c a 2 c 



dv\ 
dr). 



C = 0, 



(8d) 



(e c /r c )(dr/de) € 



B 



N c {l-al) 

N c a 2 c + (Q c /T c )(dT/dO) 
N c a c 



r c (r c - 2) 



(8f) 
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and 



C 



N c a 2 c + (Q c /T c )(dT/dQ) c 
N r r r 



-2a 2 Jl 



Equation (8d) has two roots. For radial flows, the roots are of the saddle type, where 
(dv/dr) c is real and (dM/dr) c is of opposite signs for the two roots. Here, M = v/a is 
the Mach number. Moreover, the two roots can be either of the acceleration type (A- type), 
where (dv/dr) c is of opposite signs, or of the deceleration type (D-type), where (dv/dr) c is 
negative for both roots. In the A-type, both the acceleration and wind flows accelerate at 
the sonic point. On the other hand, in the D-type, only the accretion flows accelerate, while 
the wind flows decelerate at the sonic point. 

By substituting the quantities at the sonic point, equations (7b) give the temperature 
gradient at the sonic point 



(dB\ 



Oc 

'Nr. 



2r,-3 



+ 



_r c (r c - 2) v c (l - ujj) 



dv 
dr 



(8h) 



Finally by integrating the equations of motion, we get the relativistic Bernoulli equation 
( iLightman et al. 1119751 ) 

5- (/ + 2Q K , Q) 

where £ is the Bernoulli parameter or is also known as the specific energy of flows. Since we 
assume adiabatic flows without heating and cooling, £ is a constant of motion. 



2.4. Procedure to get global solutions 



Combining equation (8b) and (5k) gives Q r in terms of r r and £. Combining it with equa- 
tion (9) gives a formula involving r c , £, and £ rtChakrabartll990l ] 1996b; F ukumura fc Kazanasl 
20071 ). If £ and £ are given, then r c is computed from the formula. Once r c is known, all the 
quantities at r c , e.g., O c , v c , (dv/dr) rc , (dQ/dr) rc , etc, are computed from equations (8a) - 
(8h). Then equations (7a) and (7b) are integrated, starting from r c , once inwards and then 
outwards, to obtain the global, transonic solutions of spherical flows around black holes. By 
this way, we can obtain two parameter {£, £) family of accretion and wind solutions. 



3. Sonic Point Properties 



In the transonic flows we study, the sonic point plays an important role. So before we 
present global solutions in the next section, we first investigate the properties of the sonic 
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point in this section. Understanding the sonic-point properties will allow us to have an idea 
of the nature of global flow structures. 

The sonic point location, r c , that is computed as a function of £ and £, is presented in 
Figure 2a. Corresponding to each set of £ and £ values, there exists a unique r c . Each curve, 
which is given as a function of £, is for a different value of £. If a flow is more energetic with 
larger £ , it is characterized by a smaller value of r c . However, at the same £ , r c is smallest 
for the e~ — e + fluid (solid line). The value of r c increases for £ < 0.2, and then starts to 
decrease for larger £. In other words, if fluids of the same £ but different £ are launched at a 
large distance away from a black hole, then the e~ — e + fluid crosses the sonic point closest 
to the event horizon, compared to the fluids of finite proton proportion. Alternatively, at 
the same r c , £ is smallest for the e~ — e + fluid, and it increases up to £ ~ 0.2 and then 
decreases for 0.2 < £ < I. Although for the same r c the e~ — p + (dotted line) is not 
most energetic, it is definitely more energetic than the e~ — e + fluid. Since fluids of different 
composition are energetically quite different at the same r c , or conversely fluids of different 
composition but the same £ form the sonic point at widely different r c , it is expected that 
the global solutions of accretion and wind flows would be quantitatively and qualitatively 
different, depending upon the composition of fluids. 

In Figures 2b and 2c, we show T c and N c as a function of r c . Equations (8b) tells that 
the sound speed at the sonic point, a c , is fixed, once r c is determined (a c implicitly depends 
on £ and £ through r c ). So plotting any variable as a function of r c is equivalent to plotting 
it as a function of a c . As noted above, at the same a c , fluids composed of lighter particles are 
colder. Therefore, in Figure 2b, at the same r c , the temperature is lowest for the — e + fluid, 
and progressively gets higher for fluids with larger proton proportions, and the maximum 
temperature is for the e~ — p + fluid. However, as noted before, higher T c does not necessarily 
ensure higher N c (i.e., more relativistic fluids). In Figure 2c, at the same r c , the e~ — e + 
fluid has the lowest N c , that is, it is least relativistic. In the range of a few < r c < 100, 
at the same r c , N c increases as the proton proportion increases for £ < 0.2, and then starts 
to decrease for 0.2 < £ < 1. This is a consequence of the competition between the thermal 
energy and the rest mass energy, as discussed in connection with Figures lc. In order to 
make the point even clearer, in Figure 2d, we show N c as a function of £ for a wide range of 
values of r c . Each curve with a single value of r c signifies fluids of different composition but 
the same sound speed at the same sonic point. N c tends to peak at some values of £, where 
the thermal contribution with respect to the rest mass energy contribution peaks. For small 
values of r c (i.e., large a c 's), a small increase of £ causes the thermal contribution to peak. 
For large values of r c (i.e., small a c 's), large proton proportions are needed to achieve the 
same. 
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As discussed in section 2.3, the roots of equation (8d) are either of the A-type or of the 
D-type. At small values of r c , the nature of the sonic point is of the A-type. It is because if 
the sonic point form closer to the central object, the flow is hotter at the sonic point (Figure 
2b), and in the wind that is thermally driven, the flow tends to accelerate at the sonic point. 
But beyond a limiting value, say r C £, the nature changes from the A-type to the D-type, 
where the wind flow decelerates at the sonic point. In Figure 3a, r c e is plotted as a function 
of £. Since at a given r c the e~ — e + fluid is thermally least relativistic, r rf is smallest for the 
fluid. The limit r rf increases with £. However, since increasing £ makes fluids 'heavy' too, 
r c e peaks around £ ~ 0.75. In Figure 3b, we plot the limiting values of £ corresponding to 
rtf, £e, as a function of £, such that for £ > En the nature of the sonic point is of the A-type, 
for £ < Si it is of the D-type. 



4. Spherical Accretion and Wind Solutions 

In this section, we present the global solutions of equations (7a - 7b) that were obtained 
with the procedure described in section 2.4. In Figure 4, we first compare typical accretion 
and wind solutions of the A and D-types for the e~ — e + fluid. The solutions of the A- 
type in the left panels have the sonic point at r c = 4, inside r c i, while those of the D-type 
in the right panels have r c = 30, beyond r C £ (Figure 3a). The accretion solutions (solid 
curves) are characterized by supersonic flows at the inner boundary and subsonic flows at 
the outer boundary (Figures 4c and 4d). The wind solutions (dotted curves), on the other 
hand, have subsonic flows at the inner boundary and supersonic flows at the outer boundary. 
The accretion flows around black holes necessarily accelerate inwards. However, the wind 
flows may accelerate (Figure 4a) or decelerate (Figure 4b) outwards. The wind solutions 
considered in this paper are thermally driven. These winds are very hot at the base, and are 
powered by the conversion of the thermal energy into the kinetic energy. It can be shown 
from equation (7b) that 



dQ 9 



2r-3 
r(r - 2) 



dv , . 

Tr <«- do) 



In other words, if the outward thermal gradient is weaker than the gravity, the outflow can 
decelerate. For the wind with r c = 30 (Figures 4b, 4d, and 4f), —d<d/dr ~ (Q/N)(2r — 
3)/[r(r — 2)] at r ~ 9.16, exactly where the outflow starts to decelerate. However, the wind 
velocity will reach an asymptotic value at r — > large, since —dQ/dr ~ (Q/N)(2r — 3)/[r(r — 
2)] ~ at large distances from the black hole. Similar relation between the gradients at the 
sonic point will determine the nature of the sonic point. It may be noted that at r c > r C £ 
(Figures 3a and 3b), such relation between (dv/dr) c and (dQ/dr) c is satisfied. Regardless 
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of accretion/wind or the type, the temperature decreases with increasing r (Figures 4e and 
4f). 

We note that the winds in our solutions are too weak to be the precursor of astrophysical 
jets, until and unless othe r accelerating processe s like those caused by magnetic fields or disc 



radiation are considered (jChattopadhyayl 120051 ) . In fact, we checked that it is not possible 



to generate the terminal speed much greater than ~ 0.8c for purely thermally driven winds, 
such as the ones that are considered in this paper. It is also to be noted that our D-type, 
wind solution is not an example of 'breeze'. A breeze is always subsonic, while the wind here 
is transonic, albeit decelerating. 

In the previous figure, we have compared the solutions with the same £ (= 0) but 
different r c . In Figure 5, we compare the solutions with the same r c (= 20) but different 
£. As shown in Figure 2a, even for the same r c , the specific energy is different for fluids of 
different £. Furthermore, the polytropic index at the sonic point is different too (N c = 1.547 
in Figure 5a, N c = 2.626 in Figure 5b, and N c = 2.271 in Figure 5c). Therefore, even if 
we fix the sonic point (and therefore a c ), the flow structure and energetics are different for 
fluids with different £. In these particular solutions, the e~ — e + fluid is not hot enough to 
drive an accelerating wind (Figure 5a), while the fluids with significant protons can do so. 
As in the previous case of decelerating wind solution (i.e., Figure 4b), in the present case the 
— e + fluid first accelerates and then starts to decelerate at r ~ 9.86. The velocity profile 
eventually tapers off to an asymptotic value at large distances away from the black hole. It 
has been shown in Figure 2b that at the same r c , adding protons increases the temperature 
at the sonic point. Larger temperature gradient causes winds of finite proton proportion to 
be accelerated at the sonic point (Figure 5b). It is seen that beyond a critical value, the 
increase in £ increases the inertia which reduces the wind speed, as is vindicated by Figures 
5b and 5c. It is to be remembered that, the D type sonic point is a reality for fluids of any 
f , provided r c > r rf . 

Although the wind solutions are noticeably different depending on £, there seems to 
be only small difference in the velocity profile of the accretion solutions. Henceforth, we 
concentrate only on accretion solutions. Such small difference in v in accretion solutions is 
expected. The accretion is generated mostly by the inward pull of the gravity, which gives 
the unique inner boundary condition for black holes, i.e., v = 1 at r = 2, regardless of other 
considerations. The pressure gradient changes the profile of v too. Since the composition 
of fluids determines the thermal state, it influences the profile of v, but the effect is not the 
dominant one. 

In Figure 6, we compare the accretion solutions with the same S (= 1.015) but different 
£. As noted below equation (9), S is a constant of motion. For r — > oo, as u t — > 1, we 
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have S — > /i^, where 



e + p 




r /+2e i 


. P . 


oo 





(11) 



is the specific enthalpy at infinity. Equation (11) tells us that at large distances from black 
holes, for the same S, T is large if £ is large. Hence, fluids with larger £ are hotter to start 
with. Therefore even for fluids with the same S, the solutions are different if £ is different. 
Figure 6a shows the velocity profile as a function of r. Here, the difference in v for fluids 
with different £ is evident, albeit not big as pointed above. Figures 6b, 6c and 6d show the 
mass density, temperature, and polytropic index. To compute the mass density, we need to 
supply the mass accretion rate, which is given as 

M = Aitr 2 u r p (12) 

from equation (6c). The mass density in Figure 6b was computed for M B = 1OM and 
M = 0.1M Edd , where M Edd is the Eddington rate of accretion. The difference in T and N 
for fluids with different £ is more pronounced. The e~ — e + fluid is slowest, densest (for 
the same M), coldest, and least relativistic. The e~ — p + fluid is more relativistic than the 
— e + fluid. But the most relativistic fluid is the one with the intermediate value of £. It 
is interesting to note that except for the e~ — e + fluid, iV is a slowly varying function of r 
for the other two fluids. Does this mean it would be sufficient to adopt the fixed T EoS with 
appropriate values of T? 

Finally in Figure 7, we compare the accretion solutions with the same temperature 
at large distances but different £. All the fluids start with T = T out = 1.3 x 10 9 K at 
r = r out = 2000. Again the mass density was computed for M B = 1OM and M = 0.1M Edd . 
It is to be noted that the fluids starting with the same T out but different £ have different 
specific energies. Hence, the velocity at the outer boundary is different too. As shown in 
Figure 7a, in these particular solutions, the e~ — e + fluid starts with a velocity substantially 
different from those of the other two fluids, so the resulting velocity profile is substantially 
different. From Figure 7d, it is clear that there are significant variations in iV for all the 
fluids. The e~ — e + fluid starts with the largest N. It is because at the same temperature, 
the e~ — e + fluid is thermally most relativistic. The behavior of N can be traced back to 
Figure la. For instance, the variations in iV tend to flatten at T > 10 10 K. In Figure 7c, for 
the fluids with £ = 0.5 and 1, T < 10 10 K for r > 100 and T > 10 10 K for r < 100. So 
significant variations are expected in iV at r>100, while the variations flatten at r < 100. 
Similar considerations will explain the variations in iV for the e~ — e + fluid. From Figure 
7d, it is clear that we need to adopt a relativistically correct EoS [equation (4c) or (4d)], 
instead of the EoS with a fixed T, in order to capture the proper thermal properties of flows 
around black holes. 



13 



In this section, we have shown that fluids with different composition can result in dra- 
matically different accretion and wind flows, even if they have the same sonic point or the 
same specific energy, or they start with the same temperature at large distances from black 
holes. So not just adopting a correct EoS, but incorporating the effects of fluid composition 
into the EoS (see equation 5e) should be also important in describing such flows. 



5. Validity of EoS 

In section 2, we have made the following assumptions for our EoS (equation 5e); fluids 
are in equilibrium, i.e., 1) the distribution of the constituent particles is relativistically 
Maxwellian and 2) the multi-components are of single temperature. However, it is not 
clear whether the conditions are satisfied. Most astrophysical fluids, unlike the terrestrial 
ones, consist of charged particles, which are collisionless, and so held together by magnetic 
fields. The constituent particles, on the other hand, exchange energies, and become relaxed 
mostly through the Coulomb interaction, which is a slow process in collisionless plasmas. In 
addition, most of the heating processes, such as viscosity and shock heating, are likely to 
affect protons. However, it is mainly the electrons which radiate. So the energy exchange 
between electrons and protons should operate, and eventually govern the thermal properties 
of fluids. 

Let t ee be the electron-electron relaxation time scale, t pp be the proton-proton relaxation 
time scale, and t ep be the electron-proton relaxation time scale. And let t pro b be the time 
scale of problem, such as the dynamical time scale, or the heating and/or cooling time scale. 
Only if t ee < ^prob an d tpp < t pro k>, electrons and protons will separately attain the Maxwellian 
distributions. And only if t ep < t pro b, electrons and protons will relax to single temperature. 

To verify the assumptions for our EoS, in this section, we compare the relaxation time 
scales with the dynamical or accretion time scale (tdyn —r/v) for an accretion solution. We 
consider the temperature range where protons are thermally non-relativistic while electrons 
are relativistic. In most our solutions in the previous section, the computed temperature 
favors this ran ge. The relativistic electron-electron interaction time scale was derived by 



Stepnev! (J1983J) 



8k 2 T 2 , . 
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The time scale for the non-relativistic proton-proton interaction is given in ISpitzerl (119621 ) 



4^ (mtfT* 
™ lnA(m p c 2 ) 3 /V T c \mj n pt ' ( ' 



-14- 



The relativistic electron-proton interaction time scale was also derived by IStepneyl (119831 ). 

1 T 



t, 



cp 



m P c A 
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We present the electron number density, n e - (Figure 8a), the three velocity v (Figure 
8b) and the temperature T (Figure 8c) of the accretion solution for the e~ — p + fluid with 
£ = 1.5247. The electron number density was computed for Mb = 1OM and M = O.lMEdd- 
In Figure 8d, various time scales are compared. All the relaxation time scales were calculated 
for the solution of single-temperature. To our surprise, it is clear that the accretion flow in the 
figure is 'too fast', such that various relaxation time scales are longer than the accretion time 
scale at least within few tens of r s . The implication of it is not clear, however. For instance, 
in relativistic plasmas, the constituent particles can be relaxed through the interactions with 
magnetic fields, too. But the relaxation will depend on the details of field configuration, 
such as the strength and the topology. Since we ignore in this study magnetic fields as well 
as other processes such as non-conservative ones, we leave this issue of the validity of our 
EoS for future studies. 



6. Discussion and Concluding Remarks 

In this paper, we have investigated the effects of fluid composition on the solutions 
of accretion and wind flows onto black holes. In order to elucidate the effects, we have 
considered a very simple model of spherical flows onto Schwarzschild black holes, and non- 
conservative processes and magnetic fields have been ignored. 

First, we have suggested an approximate EoS for multi-component fluids in equation 
(5e), and studied the thermal properties of fluids with the EoS. Three temperature ranges 
have been categorized; for kT < m e c 2 , any type of fluids are thermally non-relativistic, for 
kT > rripC 2 , any type of fluids are thermally relativistic, and for m e c 2 < kT < m p c 2 , the 
degree to which fluids are relativistic is determined by the composition of the fluids as well 
as the temperature (Figure la). Then we have shown that although at the same temperature 
the e~ — e + fluid is most relativistic (Figure la), at the same sound speed it is least relativistic 
(Figure lc), compared to the fluids with protons. It is because whether a fluid is relativistic 
or not depends on the competition between the thermal energy and the rest mass energy of 
the fluid. 

The thermal properties of fluids carry to the sonic point properties. The sound speed 
at the sonic point, a c , explicitly depends only on the sonic point location, r c (it implicitly 
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depends on the specific energy, £, and the proton proportion, £, through r c ). Therefore, 
comparing the thermodynamic quantities at the same r c is equivalent to comparing those 
quantities at the same a c . We have shown that at the same r c , the e~ — e + fluid is least 
relativistic, and a fluid with a finite £ is most relativistic (Figures 2c and 2d). 

Then, we have presented the global solutions of accretion and wind flows for the same 
r c but different £, for the same £ but different £, and for the same T at large distances from 
black holes but different £. In all the cases, the flows can be dramatically different, if the 
composition is different. This asserts that the effects of fluid composition are important in 
the solutions, and hence, incorporating them properly into the solutions through the EoS is 
important. 

Lastly, we have noted that the EoS in equation (5e) is based on the assumptions that 
the distribution of the constituent particles is relativistically Maxwellian and the multi- 
components are of single temperature. However, at the same time, we have pointed out that 
while the Coulomb relaxation times are normally shorter than the dynamical time far away 
from black holes, close to black holes they can be longer. It means that close to black holes, 
the assumptions for the EoS can be potentially invalidated. The implication of it needs to 
be understood, and we leave further consideration of this issue for future studies. 
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Fig. 1. — (a) Polytropic index as a function of the temperature, (b) sound speed as a 
function of the temperature, and (c) polytropic index as a function of the sound speed, for 
multi-component relativistic fluids of different composition with proton proportions of £ = 
(solid), 0.2 (dashed), 0.4 (long dashed), 0.6 (dashed-dotted), 0.8 (long dashed-dotted) , and 
1 (dotted). Hereafter, the temperature in figures is given in units of Kelvin. 
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Fig. 2. — (a) Sonic point location as a function of the specific energy, (b) temperature at the 
sonic point as a function of the sonic point location, (c) polytropic index at the sonic point 
as a function of the sonic point location, for transonic flows of fluids with £ = (solid), 0.2 
(dashed), 0.4 (long dashed), 0.6 (dashed-dotted), 0.8 (long dashed-dotted) , and 1 (dotted), 
(d) Polytropic index at the sonic point as a function of proton proportion for the flows with 
the sonic point location of r c = 5 (solid), 105 (dotted), 205 (dashed) and 13505 (long dashed). 




Fig. 3. — Limiting values of (a) the sonic point location and (b) the specific energy, which 
divide the domain of the A-type sonic point roots from that of the D-type sonic point roots, 
as a function of proton proportion. 
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Fig. 4. — Examples of transonic accretion (solid) and wind (dotted) solutions of the A-type 
with the sonic point at r c = 4 (£ = 1.3, N c = 2.163) (left panels) and of the D-type with the 
sonic point at r c = 30 (£ = 1.0016, N c = 1.548) (right panels). The e~ — e + fluid (£ = 0) 
is considered. The radial three- velocity (top), Mach number (middle), and temperature 
(bottom) are shown as a function of radius. 
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Fig. 5. — Comparison of accretion (solid) and wind (dotted) solutions with the same sonic 
point location of r c = 20 but different proton proportions of (a) £ = (£ = 1.0039), (b) 
£ = 0.5 (£ = 1.0337), and (c) £ = 1 (£ — 1.0239). The radial three-velocity is shown as a 
function of radius. 



-22 - 





1 2 3 1 2 3 

log(r) log(r) 



Fig. 6. — Comparison of accretion solutions with the same specific energy of S = 1.015 but 
different proton proportions of £ = (solid), £ = 0.5 (dashed), and £ = 1 (dotted). The 
sonic point locations are r c = 11.0 for £ = 0, r c = 38.267 for £ = 0.5, and r c = 28.972 for 
£ = 1. The radial three- velocity (a), mass density (b), temperature (c), and polytropic index 
(d) are shown as a function of radius. The mass density was computed assuming the black 
hole of Mb = 1OM and the accretion rate of M — 0.1MEdd 5 and is given in units of g cm -3 . 
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Fig. 7. — Comparison of accretion solutions with the same temperature of T out = 1.3 xlO 9 K 
at the outer boundary r out = 2000 but different proton proportions of £ = (solid), £ = 0.5 
(dashed), and £ = 1 (dotted). The sonic point location and the specific energy are r c = 3.5 
and 8 = 1.6322 for £ = 0, r c = 333.3 and 8 = 1.0008 for £ = 0.5, and r c = 806.4 and 
8 = 1.0002 for £ = 1. The radial three- velocity (a), mass density (b), temperature (c), and 
polytropic index (d) are shown as a function of radius. The mass density was computed 
assuming the black hole of Mb = 1OM and the accretion rate of M = O.lMEdd, and is given 
in units of g cm' 3 . 
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Fig. 8. — (a) Electron number density, (b) radial three-velocity, and (c) temperature as 
a function of radius in an accreting flow of the e~ — p + fluid (£ = 1) with £ = 1.5247. 
M = O.lMEdd onto a black hole of Mb = 1OM was assumed, (d) Comparison of various 
time scales of the same flow. Different curves represent the accretion time td jn (solid), the 
electron-electron relaxation time t ee (dotted), the proton-proton relaxation time t pp (dashed), 
and the electron-proton relaxation time t ep (long-dashed). 



